Effects of an external magnetic field on the gaps and quantum corrections in an 
ordered Heisenberg antiferromagnet with Dzyaloshinskii-Moriya anisotropy 
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We study the effects of external magnetic field on the properties of an ordered Heisenberg antifer- 
romagnet with the Dzyaloshinskii-Moriya (DM) interaction. Using the spin-wave theory quantum 
correction to the energy, on-site magnetization, and uniform magnetization are calculated as a func- 
tion of the field H and the DM anisotropy constant D. It is shown that the spin-wave excitations 
exhibit an unusual field-evolution of the gaps. This leads to various non-analytic dependencies of the 
quantum corrections on H and D. It is also demonstrated that, quite generally, the DM interaction 
suppresses quantum fluctuations, thus driving the system to a more classical ground state. Most 
of the discussion is devoted to the spin-S, two-dimensional square lattice antiferromagnet, whose 
5=1 case is closely realized in K2V3O8 where at H = the DM anisotropy is hidden by the 
easy-axis anisotropy but is revealed in a finite field. The theoretical results for the fleld-dependence 
of the spin-excitation gaps in this material are presented and the implications for other systems are 
discussed. 

PACS numbers: 75.10.Jm, 75.30.Ds, 75.40.Cx, 75.40.Gb 



I. INTRODUCTION 

The studies of quantum antiferromagnets (AFs) in ex- 
ternal magnetic field attract significant attention because 
these systems exhibit a variety of unusual quantum- 
mechanical phenomena of general interest. The modern 
high magnetic field technology, materials' synthesis, and 
experimental probes allow precise measurements in the 
regime of the field strength comparable to the charac- 
teristic exchange constant of a system. This has made 
possible the observations of condensation of triplet exci- 
tations in a variety of chain, ladder, and weakly-coupled 
dimer compounds p^^ii magnetization plateaux in frus- 
trated magnetS)^ and other new effectsi^ It turned out 
that in many cases experimental data deviate signifi- 
cantly from the theoretical predictions based on the pure 
isotropic Heisenberg model in external field&L^ Such 
deviations are due to anisotropics, most notably the 
Dzyaloshinskii-Moriya (DM) anisotropy, which are usu- 
ally small and often neglected from zero-field considera- 
tions. Not only such anisotropics can induce qualitatively 
new effects, but also the strength of such effects seems to 
be significantly amplified by an external field. In par- 
ticular, in the one-dimensional (ID) spin-i chains with 
the DM interaction spin excitation spectrum exhibits a 
gap A oc {DH)^/^ in contrast with the pure Heisenberg 
case where the spectrum remains gapless^*^ Thus, while 
usually treated as an insignificant contamination of the 
problem, the DM interaction can give rise to an interest- 
ing phenomena on its own. 

In the case of the ID spin-chains an unusual field- 
dependence of the gap has been first documented 
experimentally^ and has received a theoretical expla- 
nation soon after tYvaXw^ For the higher-dimensional 
cases there is a recent theoretical expectation^'"'^^'^^ that 
the field-dependence of one of the gaps should be A cx 
{DHy/^. However, such a behavior has not been ob- 



served in any known system yet. There is also a dis- 
agreement between the recent theoretical worksiiii^ in 
the low-field regime with the earlier low-field studies of 
the Heisenberg-|-DM problem li^ii^ii^ The latter works 
predict a non-zero gap and also weak ferromagnetismi^ 
in zero field, while the former do not. We will show below 
that this discrepancy comes from the higher-order con- 
tributions, which yet result in the effects of order 0{D) 
in the excitation gap and spin canting angle. We will 
also show that in a realistic system the low-field behav- 
ior can be more delicate. The field-dependence of the 
gap, quite generally, takes the form A oc \J D{H — H*), 
where H* ~ D can be both positive and negative depend- 
ing on the system. This clarifies theoretical expectation 
for experiments. 

Furthermore, most of the studies have excluded from 
consideration the quantum effects associated with the 
DM interaction. The generic question is: does the DM 
interaction enhance or suppress quantum fluctuations? 
This is important for the systems close to the spin-liquid 
state, such as triangular lattice AF CS2CUCI4 where the 
DM term is significant ji^iiSiiS and also for the frustrated 
spin systems where quantum effects may help to select 
the ground statei^SiSi 

In this work we study the excitation spectrum and var- 
ious T — static properties of the square lattice, spin-S* 
Heisenberg model with the DM interaction. This is the 
simplest case of a higher-dimensional system in which 
generic trends of the mutual effect of the DM interaction, 
external field, and quantum fiuctuations on the proper- 
ties of an ordered AF can be studied. This choice is also 
motivated, in part, by the recent experiments in K2V3O8, 
a 2D, S — ^ AF, which demonstrates an unusual spin- 
reorientation transition in a small fieldiS^ This effect has 
been understood as resulting from an interplay of the DM 
and easy-axis anisotropics. 

Altogether, the purpose of this work is threefold. First, 
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to reconcile the earlier low-field and recent high-field the- 
oretical studies of the Heisenberg-|-DM-|-field problem. 
Second, to study the dependence of the quantum cor- 
rections on the DM interaction. Third, to analyze an 
experimental system to which these results can be ap- 
plied. 

We find that the gap in one of the magnon branches 
shows the following behavior when H _L D: A^ttt — D 
in H ^ field, A^^ = VD for H ~ 0.7H„ and 
A^TT — V^D"^/^ for H = Hs, where Hs is the satura- 
tion field and the energy units = AS J are used. These 
findings are in agreement with the other works >ii*iSii^ 
We also show that, for the fields < H < Hg, the 
DM interaction leads to the suppression of the quantum 
fiuctuations. The dependence of the quantum correc- 
tions on D for various quantities is often non-analytic. 
It is only in the regime of the field close to the satu- 
ration, H ~ Hs, where the DM term enhances quan- 
tum fluctuations, effectively leading to a proliferation of 
the quantum effects into the classical saturated phase 
{H > Hg). We also show that some quantum correc- 
tions are singular in the limit of both H and D going 
to zero, namely lim//^o limri^o 7^ liniu^o liniff-^o- We 
demonstrate that K2V3O8 is an excellent candidate for 
the observation of the unusual field-dependence of the 
gap, characteristic to the 2D and 3D AFs. In fact, an 
additional easy-axis anisotropy helps to switch between 
the Ising-like and easy-plane behavior in this system, 
making it potentially possible to observe a non-analytic 
A (X yjH — H* behavior of the gap.^^ 

The paper is organized as follows. In Sec. ^Iwe de- 
scribe the technical details, derive the spin-excitation 
spectrum, and present the results for the quantum cor- 
rections to various quantities in different field regimes. 
In Sec. mil we give a quantitative discussion of the ex- 
citation spectrum in K2V3O8. We conclude by Sec. IIVI 
which contains a brief discussion of our results. 



II. MODEL, SPECTRUM, AND QUANTUM 
CORRECTIONS 

As is mentioned in Introduction, we would like to con- 
sider the simplest possible case of an ordered AF with 
the DM anisotropy in an external field, which would al- 
low to investigate the role of quantum corrections in it. 
We therefore study the spin-S", nearest-neighbor Heisen- 
berg model on a square lattice with the DM interaction: 




{ij) denotes summation over bonds and H is the mag- 
netic field in the energy units {gus = !)■ The direc- 
tion of the (staggered) DM vector will be chosen along 
the z-axis, D^^ = (— 1)^(0, 0, D). Aside from being the 
simplest and the most commo n^ ^1^^ choice, this Hamil- 
tonian also corresponds closely to K2V308>^ We would 



like to note that, generally, this form of the Hamiltonian 
is incomplete, because, as noted in Refs. ■25..26...27„ the 
same microscopic processes that generate the DM cross- 
product term also yield the higher-order 0{D^) terms, 
which contribute to some observables on equal footing. 
Without going into much details we refer to the same 
.^pj.jj.g25j26j27 -^^jjich show that these additional terms can 
be cast in the form of an easy-axis anisotropy with the 
"easy" axis along the D-vector, also the case realized 
in K2V3O8. Quite generally, the DM constant D and 
the easy-axis anisotropy constant E are not related to 
each other in a straightforward way and thus can be 
treated as phenomenological parameters. Therefore, in 
the following we would like to concentrate on the most 
commonly used Hamiltonian which only contains the 
cross-product DM-term and will consider the effect of 
the easy-axis anisotropy later in the context of K2V3O8, 
Sec. mil It is convenient to use units J = 1 in Eq. 
and this convention is used throughout the paper unless 
noted otherwise. 



A. Classical limit 

It is important to discuss the classical limit of the 
model ^ first, as for the consideration of the excitation 
spectrum the quantization axes for spins will be chosen 
along their classical directions. The effect of the DM in- 
teraction is twofold: it makes energetically favorable for 
the spins to stay in the plane perpendicular to the di- 
rection of D and, without the J-term, it would make the 
spins in different sublattices to be under an a = 7r/2 angle 
to each other (in the notations of Fig. ^that corresponds 
to the angle = (tt — a)/2 = 7r/4). While the former 
trend can only be hindered by additional anisotropics 
(see next Section) the latter always competes with the 
much larger Heisenberg term which prefers the spins to 
stay antiparallel. Thus, at _ff = in the ground state of 
the model Q spins lie in the x — y plane and are canted 
towards each other. According to Fig. ^this configura- 
tion has a finite angle ip{H = 0) = i tan~^ D resulting in 
a non-zero uniform magnetization M = SsmLp ~ SD/2 
(per spin) . Therefore, the H — limit of the model ^ 
exhibits an effective easy-plane anisotropy and a weak 
ferromagnetism, a phenomenon the DM interaction was 
originally proposed for^ 

It is clear now that there are two distinct choices for 
the direction of the external field: one is along the DM 
vector, H |j z, and the other one is perpendicular to it, 
H _L z. We will be mostly concerned with the latter case 
of the in-plane field as it leads to most non-trivial results 
and will only briefly discuss the former. From the sym- 
metry point of view these two cases are clearly distinct 
because the symmetry is lowered by the DM interaction 
already in zero field from the full rotational 0(3) sym- 
metry to the easy-plane 0(2) symmetry. Then the field 
along the z-axis will only cant the spins towards itself 
without affecting the freedom of choice of the sponta- 




FIG. 1: (Color online) (a) H||D||2: configuration. The canting 
angle f in the x—y plane is induced by the DM interaction and 
does not depend on the field. The canting angle d towards the 
field (z-axis) is given by Eq. J^. (b) H _L D configuration. 
The canting angle is given by Eq. Q. 



neous in-planc spin direction, thus leaving 0(2) intact. 
On the other hand, the in-plane field will break the re- 
maining continuous 0(2) symmetry and "pin" the spins' 
direction in the plane. 

Canting angles at D = 0. Since the magnitude of the 
DM interaction is generally small one may wonder how 
much difference its presence makes for the ground state 
of the model (Q) for the finite, and not necessarily small, 
values of H. While we will address this issue in detail 
below it is instructive to consider the spin configuration 
when D = 0. If we neglect the DM interaction in Q 
the only effect of the external magnetic field is to induce 
a finite uniform magnetization by canting spins towards 
the field. The canting angle is i/jq = sm^^{H/8S) and 
the spins become fully polarized in the saturation field 
Hs = 8S (where 8 = 2z, z = 4 is the coordination number 
for the square lattice). 



1. Out-of-plane field. 

We now consider the field-dependence of the canting 
angles in the presence of the DM interaction. If the field 
is perpendicular to the x — y plane the canting induced 
by the _D-term and the one induced by the field lie in the 
orthogonal planes, see Fig. Q^a). Therefore, one should 
expect the DM-induced in-plane canting angle tp to be 
independent of the field. Using the angle notations of 
Fig. n^a) the classical energy of the model Eq. can 
be written as: 



E,. 



1 - sin^ 9 (2 cos^ + D sin 2(p) 



H 
2S 



cos 9 , (2) 



minimizing which gives the canting angles, N is the num- 
ber of lattice sites. As expected, f — ^ tan~^ D is inde- 
pendent of the field and the field-induced out-of-plane 
canting angle should be found from: 



H 



cosf = 



85* 1 + VTT^ 



(3) 



Thus, the only effect of the small DM interaction for the 
case of H||D is the slight change of the saturation field 
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FIG. 2: (Color online) The difference of the in-plane canting 
angle <^ and its _D = limit ifo v.s. field in the H _L D con- 
figuration for two representative values of D = 0.001 (dashed 
lines) and D = 0.04 (solid lines), D is in units of J. Inset: 
the canting angle v.s. field. The asymptotic behavior of 
(/5 — (/5o = —DHs/2{H — Ha) is shown by the dotted line. 
D^^^ values for both choices of D are shown by the horizontal 
dashed lines. 



from = 85 to Hi = 45 (l + ViTD^) . 



2. In-plane field. 



For the in-plane field the changes for the canting angle 
are less trivial. The spins are in the x — y plane, see Fig. 
^b). From minimizing the classical energy: 



— COS 2ip — D sin 2ip sin ip , 



2NS^ 

the canting angle obeys the following equation 

D cos2lp H 
2 cos (p 8S 



sm if ■ 



(4) 



(5) 



Without the DM interaction Eq. jSJ yields a familiar 
i/-dependence of the canting angle: (po = sm^^{H/8S) 
(0 in the subscript refers to D = 0). To demonstrate 
the effect of the DM interaction on the canting angle in 
the H _L D case we plot the difference of {ip — ipo) v.s. 
field for two representative values of D in Fig. |21 One 
can see that the difference is D/2 for small field and it 
vanishes for the field H ~ 0.707Hs. This field makes 
the angle between the spins to be 7r/4 for any value of 
D. Thus, the DM-term in the energy is always fully 
minimized for such a field. The most striking feature in 
this plot is the sharp singularity at _ff = H^. This is 
due to the fact that the D — canting angle (po exhibits 
a kink at Hs where it reaches the maximal value tt/2. 
In contrast, the field-dependence of the canting angle ip> 
is smooth through the saturation field, as shown in the 
inset of Fig. |21 In fact, the spins never fully saturate, 
only as H oo. The asymptotic behavior of the angle is: 
(p-n/2 = -DHs/2{H-Hs), which is shown in Fig. Elby 
the dotted line. Therefore, the saturation field formally 



ceases to exist for the H _L D case. In the rest of the 
paper we refer to the "bare" {D = 0) value Hg = 8S as 
to the value of the saturation field. It is interesting to 
note that the magnitude of the deviation of the canting 
angle from its D = value is quite significant close to the 
saturation field. Taking H = Hg and assuming D <^ 1 
one can find from Eq. that cosip ~ it /2 — ~ D^l^ , 
a value which is not necessarily small even if D itself is 
reasonably small. The D^l^ values for the choices of D 
are shown in Fig. |21by the horizontal dashed lines. 

B. Spectrum and gaps 

Technical approach. To apply the spin- wave theory to 
the model (^3) we follow the general prescription of Ref. 
I23 . Namely, the rotating local reference frame for each 
sublattice is introduced in which the quantization axis zq 
is directed along the classical spin direction. In such a 
local frame the standard bosonic representation for the 
spin operators is used: 

= S- a^a, 

5- = ^^^i^ = /f(a + at), (6) 

where we omit higher-order boson terms from thus 
restricting ourselves to the linear spin-wave approxima- 
tion. Higher-order, as usual, means both the higher num- 
ber of boson operators and the higher powers of 1/S. 
Then one needs a transformation of the spin components 
from the local reference frame to the laboratory reference 
frame. Such transformations are again based on the clas- 
sical considerations of the canting angles similar to the 
ones given in Sec. Ill Al and they depend on the relative 
direction of the field, anisotropics, etc. To give a specific 
example we provide here such a transformation for the 
H _L D configuration for the vectors D||z and H||t/, see 
Fig. Hb): 

= -5f ° sin + i-iy cos if 
Sf = 5f sin + (- 1)* cos (7) 

QZ qVO 

Using relations of that type and representation © one 
can rewrite the Hamiltonian in the form: 

Ti = Ed + T-Llswt + . . . , (8) 

where Ed is the classical energy discussed above, the 
ellipsis stands for the higher-order terms, and the spin- 
wave Hamiltonian can be written in the general form: 

Ulswt = 45 ^ (Ci + C27k) ala^ (9) 
+ (4a-k + "kfl-k) ^ , 




k/ll in (1,1) direction 

FIG. 3: (Color online) Linear spin-wave result, Eq. 1131 . 
for the magnon spectrum in a field (D = 0), for H = 0, 
H = Hs/V^, and H = Hs along the main diagonal of the BZ. 
Inset: full BZ with the direction of the cut. Magnetic BZ is 
shown by the dashed diamond. 

where 7k = (cosfc^: + cosfcj,)/2 and the summation over 
k is over the full Brillouin zone (BZ). All the coefficients 
Ci are the functions of the field and anisotropics. As is 
noted in Ref. (2^ the minimization of the classical energy 
ensures the absence of the terms linear in a, or in the 
expansion of the Hamiltonian Then, the spin- wave 
Hamiltonian is readily diagonalized by the Bogolyubov 
transformation and takes the form: 

nLSWT ^SE + J2 Wk^ftk ' (10) 

k 

where the magnon energy is 




and 

'^^= ^E(^k-45Ci) . (12) 

k 

As an example of the use of such a procedure, one can 
neglect the DM-term for a moment and, after some alge- 
bra, obtain: Ci = 1, C2.3 = (— cos 2(/jo ± l)/2. This leads 
to the following expression for the spin-wave spectrum: 

c^k = 45V(l + 7k)(l-cos2(^o7k) , (13) 

where cos2(po = 1 ^ "^(H/Hs)^, in agreement with the 
results of Ref. HI 

General features of the spectrum. Before we discuss 
the details of the gaps' dependence on both the field and 
the DM interaction, we would like to outline the broad 
features of the spin-wave spectrum. Our Fig. |3| shows 
the D = spectrum given by Eq. H13|) along the main 
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FIG. 4: (Color online) Linear spin-wave result, Eq. I|14|l . for 
the magnon spectrum in zero field for D = 0.04, D — 0.4, and 
D = 1.0 along the main diagonal of the BZ. 



diagonal of the Brillouin zone (see inset) for three rep- 
resentative values of the field. In the H = case the 
spectrum has two gapless modes, at (0,0) and (tt, tt). In 
the field one of the modes develops a gap, Aqo = H, 
which is strictly equal to the value of the field. This 
mode corresponds to the uniform precession of the field- 
induced magnetization about the field directioniS The 
other mode must remain gapless as it corresponds to the 
Goldstone mode related to the spontaneous breaking of 
the remaining 0{2) symmetry of the spins in the plane 
perpendicular to the field. 

Two remarks are in order. First, while it is conve- 
nient to use the "full" BZ language here with a single 
magnon branch, one can employ an equivalent picture 
of two magnon branches within the magnetic BZ (shown 
in the inset of Fig. 13 by the dashed diamond) . The 
neutron scattering will observe these branches near the 
(tt, 7r)-point in the different components of the structure 
factor: 5^ oc 6{u} — Wk) and S^^ oc 5{u! — Wk-(7r,7r))j 
see Ref. Second, we neglect the higher-order non- 

linear terms in the spin-wave Hamiltonian. While they 
do not change most of the results qualitatively, in the 
fields close to the saturation they induce an instability 
of the single-magnon branch towards decays into a two- 
magnon continuum. Since such an instability does not 
concern the k = (0, 0) and k = (tt, tt) modes we do not 
take this effect into account in this work:-° 

Another general feature of the spin- wave spectrum be- 
havior can be studied by considering the case H = 
and D ^ Q. Using Eqs. ©, {7)), and ® one obtains 
Ci = VTTD^, C2,3 = {-Vl + D^ ± l)/2, which gives: 



;^="=45v/rT:D^W(l-7k) 1 



7k 



,(14) 



This spectrum is plotted in Fig. ^ for three different 



values of D. Again, one of the modes is gapless due to 
the 0(2) symmetry, which is now related to the spon- 
taneous choice of the direction of the DM-induced fer- 
romagnetic magnetization vector. The gap in the other 
mode is Atttt — 4:SD for small D. These results are very 
similar to the ones obtained in Ref. in the context of 
La2Cu04. 



1. Out-of-plane field. 

Having outlined the general trends of the spin-wave 
spectrum behavior we return to the original problem of 
the DM interaction in a field. We first consider H||D||z 
configuration. Recalling our discussion about the canting 
angles for this case in Sec. Ill Al 1 . and using Eq. (PJ, after 
some algebra one obtains: 



(15) 



C2 — Ct. 



-1 



l + \/l+D2 

This gives the following spectrum: 



ul = AS^l + iy^ V(l+7k)(l-7k'5) , (16) 

with.^ri- ^ ' 



i + VTTd^J VTTd^' 

As is discussed in Sec. Ill Al this configuration of the field 
and the DM vector preserves the 0(2) symmetry of the 
problem. Because of that the spectrum remains essen- 
tially identical to the field-only, D = case. Fig. Eland 
Eq. (|13|l . That is, the (tt, 7r)-mode is gapless and the 
(0, 0)-mode is gapped with Aqo = ASy/D^ + A{H/H,)^. 
The other changes concern the saturation field, which 

is increased to = 45* + + D"^), and the total 
width of the spectrum, which is increased by the factor 
Vl + D'^. With these results we conclude our considera- 
tion of the HI ID case. 



2. In-plane field. 

We finally arrive to the most interesting part of the 
consideration: the orthogonal configuration of the field 
and the DM vector, H _L Djjz. As we have seen in Sec. 
Ill AI 2. this case is characterized by the absence of the 
true saturation of the spins and by an intriguing behav- 
ior of the canting angle. As is already mentioned, the 
non-zero DM constant and external field break all con- 
tinuous symmetries and thus no excitation branch should 
remain gapless. However, since the symmetry-breaking 
effect due to a small DM interaction should be weak one 
expects the gap associated with it to be small as well. 

With these ideas in mind one can use Fig. ^b), expres- 
sion for the canting angle for this case, Eq. jSJl, transfor- 
mation from the local to the laboratory frame, Eq. 0, 



and the transformation of Eq. Q to Eq. (|10|l to obtain: 

(17) 



Ci = 1 + D tan if , 

C2 + C3 — ~ cos 2ip — D sin 2ip , 



where the canting angle is defined from Eq. (O. This 
leads to the following expression for the spin-wave dis- 
persion: 

cot = 45^ (Ci + 7k) (c^ - - i| sin^^ .(18) 

Just to verify the consistency of this result with our pre- 
vious considerations one can set D = and obtain from 
Eq. ©: Ci = 1 and sin if = H/8S. This yields the 
result for the Heisenberg+field spectrum obtained in Eq. 
(|13|l and shown in Fig. |2| On the other hand, setting 
field to zero yields tan2i^ = D from Eq. ^ and, after 
some algebra, one reproduces the H = spectrum for the 
Heisenberg-|-DM problem given in Eq. H14|) and shown 
in Fig. H 

{Q,0)-gap. Given the above discussion, the broad 
features of the spectrum uj^ should look like a convolu- 
tion of the Figures 13 and 21 Namely, it should have gaps 
at both (0, 0) and (tt, tt) points that evolve with field. We 
first study such an evolution for the gap at k = (0,0), 
which is given by: 



Ann = 45", 



25' 



sinifi ( 1 + —Dta.n(p 



(19) 



with H- and D-dependence of ip to be defined from Eq. 
©. This gap is plotted as a function of field in Fig. 
|Slfor two representative values of D. We recall that in 
the D = case Ago grows linearly with H with a slope 
equal to unity. Since we plot Aoo/4S' v.s. H/Hg with 
Hs = 85, this translates into the slope = 2 in these units. 
One can see that for the small values of D the difference 
of our results from D = Q case mainly concerns small 
fields H <^ Hg- In this limit, assuming I? ^ 1, one can 
simplify Eqs. H19(l and (0) to obtain: 



A, 



00 



85. 



H f H 



1 



-D 



Hs \Hs 2 
y/ISD for < L> , 

H + 2SD for D<^H <$:Hs 



(20) 



Thus, for very small fields this gap is (x V H and for 
H > D it is a linear function of H with an offset oc D. 
This behavior can be seen in the inset of Fig. [S] which 
zooms into the region of small fields. Altogether, the 
changes from the DM-term in the field-dependence of the 
k — (0, 0) gap are small and are due to the fact that 
the DM interaction induces a non-zero canting of spins 
already in zero field. 




H/H^ 

FIG. 5: (Color online) Aoo/4SJ v.s. H/Hs for D = 0.001 
and D = 0.04. Inset: same plot in the region of small fields. 

{Tr,Tr)-gap. The evolution of the k — (tt, tt) gap is 
much more interesting, because it corresponds to a mode 
which would be gapless in D = case. Using Eqs. (|18f) 
one arrives to: 



A^y, ^ 45W2Dtan(^ 1 



Dtantp- g^sin<y5 ) , (21) 



which clearly vanishes as D ^ 0. The angle ip is to 
be defined from Eq. Q as before. Considering D ^ 
1 limit one can identify three different regimes for the 
field: H ^ D, Hs > H ::i> D, and H ~ Hs. After 
some algebra one can simplify Eq. H21|) and obtain the 
following expressions for the gap in these regimes: 



A, 



^/ASD ^/H + ASD 



VISD 



\ 



H^ 



H 

85 



45\/3i?2/3 



for H ^ D 



for H:$> D , (22) 



for H = H, 



These trends are clearly seen in our Fig. for a small 
value of D = 0.001. The gap evolves from A^7r/45 = D 
in zero field to its maximal value A^7r/45 = D^^^ at 
H ~ 0. 707 Hs, and then has a minimum at Hg where it 
reaches the value \/3D^^^ (marked by the dashed hori- 
zontal lines in Fig. The origin of the maximum can be 
simply understood from Eq. (|21|) assuming H ^ D and 
neglecting D from everywhere except the common pref- 
actor. This gives another expression for the (tt, 7r)-gap in 
the intermediate-field regime: 



Atttt — 45 D sin 2(p, with sin 1^9 w — 

85 



(23) 



This expression obviously has a maximum at ip = 7r/4 
{H — Hs/\/2), the canting angle at which the DM inter- 
action is fully satisfied. 
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H/H^ 

FIG. 6: (Color online) Field dependence of Ann/4:SJ v.s. 
H/Hs for D = 0.001, Eq. lUTJ. The gap values at H = 0, 
H — Hs/V^, and H = Hs are marked by the horizontal 
dashed lines. Dotted lines are: for H < Hs - the low- and 
intermediate-field approximations for the gap, Eq. 122L for 
H > Ha - the D = gap dependence Atttt = H — Ha. Arrows 
mark H = Ha/y/2. Inset: same for D = 0.04. 



One can see that for the higher value of D = 0.04 
the maximum and minimum in the field-dependence of 
the gap are only weakly pronounced (inset of Fig. IBJ, 
while the zero-field gap = D persists. Note that the zero- 
field gap is a salient feature of the model (Q), intimately 
related to the phenomenon of the weak ferromagnetism. 
The origin of the minimum and of the unusual 2/3 power 
of D can be traced to the classical consideration we gave 
for the canting angle at the saturation field in Se. Ill AI 2. 
At H — Hs we have cos ip ~ D^^^, which yields Atttt oz 
£)2/3 fj-Qjjj H21|l . For an interesting discussion of that 
fractional power see Ref. [T^i 

The minimum in AT^-rr{H) disappears, or, more pre- 
cisely, merges with the maximum and becomes an inflec- 
tion point, at a somewhat higher value of — 0.066. 
A similarly small number has been obtained in Ref. |2 
in the mean-field study of the Lieb-Mattis modeli^ It is 
surprising to have such a small "magic number" and one 
may wonder whether there is some hidden small param- 
eter behind it. It turns out that it is a small difference 
of the powers of D in A™^^ and A™" which is respon- 
sible for that smallness. One can estimate the "criti- 
cal" value of D at which the minimum disappears as the 
point where A™''^ « A™". This gives: dI^^ = 1/V3, 
or Dc — 1/27 ~ 0.037. The actual number is about two 
times higher, but already at JD = 0.04 the minimum has 
almost faded away, see Fig. El 

To conclude the description of Fig. ^ the dotted 
line for H < Hg is the convolution of the low- and 
intermediate-field results from Eq. (|22|l . which gives a 
very close description for the gap behavior up to H ^ Hg . 
The dotted line for H > Hgis the D = field dependence 



of the gap above the saturation field: Aj^t^ = H—Hg. One 
can see that the gap merges with this asymptote quickly. 
We add that the results similar to that in Fig. Elwere 
also reported for a ID model with the DM interaction in 
Ref. [lE 

Our result for Ago and A^ttt in the low-field regime, 
Eqs. H20I) and (|22l) . arc in agreement with the results of 
the earlier works from the 1960sii^ For the intermediate 
and high-field regimes our results in Eq. (|22|l agree with 
the results of the recent works on the coupled chainsii 
and Lieb-Mattis modelii^ These works used the effective 
staggered field in place of the DM interaction that allows 
to simplify the model but introduces an approximation 
that neglects effects of order 0{D) in the gaps. Such an 
approximation has created a discrepancy of these works 
with the earlier studies in the low-field regime. Namely, 
these recent works predict no weak ferromagnetism and 
no gap for the (tt, 7r)-mode in zero field. Our results for 
Atttt, Eqs. (|21|l and provide a reconciliation of the 
earlier low-field and recent high-field studies. We also 
have provided a detailed analysis of the gap in the other 
spin-wave branch, the (0, 0)-mode, Eqs. H19|l and H20|l . 
which is gapless in zero field and has a oc behavior in 
a small field regime that can be mistaken with the results 
of Refs. [miil for the (tt, 7r)-mode. 

This unified description of the behavior of the gaps 
in external field is also important as it clarifies the ex- 
perimental prediction for the gaps' dependence on the 
field. In the systems with the DM anisotropy which 
can be described by Eq. instead of the literal 

cx \/ DH non-analytic behavior of the (tt, 7r)-gap one 
should expect a non- analytic increase from a finite gap, 
cx ^JD{H + ASD), according to Eq. As we will 

show in Sec. IIIII in real systems this behavior is further 
altered by the presence of other anisotropics. 

C. Quantum corrections. 

We gave an extensive discussion of the (tt, 7r)-gap in 
Sec. Ill Bl because it will be also responsible for many of 
the changes in the behavior of quantum corrections due 
to the DM interaction. 

As is mentioned in Introduction, very little system- 
atic discussion of the influence of the DM interaction 
on the quantum effects in the Heisenberg-like systems 
has been given in the literature. The recent exception, 
Ref. [T3. considers the triangular-lattice AF-|-DM-|-field, 
to model the existing experiments in CS2CUCI4. How- 
ever, the results for that system are severely complicated 
by the other effects such as the field-dependence of the 
AF ordering vector and are hard to generalize to the other 
systems. Considering a much simpler, non-frustrated AF, 
we would like to ask a simple question: whether the quan- 
tum fluctuations are enhanced or suppressed by the DM 
interaction. 

Formulae. We are going to answer this question by 
analyzing the field-dependencies of the 1 / S'-corrections 




FIG. 7; (Color online) l/5'-corrections to the ground-state 
energy (left axis) and the on-site magnetization (right axis) 
for L» = (dashed) and D = 0.04 (solid) v.s. H/Hs, S = 1/2. 



to the ground-state energy, 5E, local (on-site) magne- 
tization, 5S, and the uniform magnetization 5M. As is 
discussed in Sec. Ill BI the out-of-plane field case, H||D||z, 
is almost identical to the D = case, considered in Ref. 
12^ Because of that, we are going to consider only the 
case of the in-plane field, H _L D||z. 

The l/S'-correction to the ground-state energy (per 
spin) is given in Eq. lO, which we rewrite here 



6E = 2Sj2('^k-Ci) , 



(24) 



using the dimensionless spin- wave energy i^k = uj^/4:S, 
Eq. H18|) , normalized to 45* to keep the S'-proportionality 
explicit. The definition of the l/S'-correction to the lo- 
cal (on-site) magnetization SS follows straightforwardly 
from: 



ss = 



= S- (atfl) =S-SS , 
/ t \ fCi + C27k 



2 ^ 

k 



- 1 



(25) 



where the field- and DM-dependent constants Ci are 
listed in Eq. ifTTjl . The classical part of the per-spin uni- 
form magnetization is simply given by the canting angle: 
Mci — Ssinip, see Fig. ^^b). A more rigorous definition: 



M 



dE 
dH 



(26) 



where E is the energy per spin from Eq. leads to 
the same answer for the classical part and also yields the 
l/S'-correction SM as follows: 



M = M, 



cl 



SM = 2Slp' 



SM ^ Ssintp- SM 

dCi 

dip 



(27) 



D = 


0.04 


— - D = 






CO 




FIG. 8: (Color online) Uniform magnetization (classical part 
-I- 1/S-correction), right axis, and the 1/5 part to the uniform 
magnetization normalized to its classical part, left axis, v.s. 
H/Hs. D = (dashed) and D = 0.04 (solid). Arrows at the 
left axis show lim//^o limn^o, Eq. I)29|l . and lim_D^o lim_ff^o, 
Eq. lOni- S = 1/2. 



where 



and 



dip 

'dH " 

dCi 

dip 



1 

8S 



cos p + D sin p) 1 



D 



2 cos^ ip 



cos^ p 



(28) 



in which we used Eq. jSjl and an explicit expression for 
Cl = 1 + D tan p from Eq. ATI) . Note that the order of 
SE is 0{S) while SS and SM are 0(1). 

General features. Having at hands expressions for 
the 1/S' corrections to E, (S^), and M, Eqs. it^ . 
and H27|) . together with the definitions of the spin- wave 
energy, Eq. (|18l) , constants Ci , Eq. H17|) , and the canting 
angle, Eq. jSJ, one can calculate their field-dependence. 
Our Figure [7| shows the results for SE/iS and SS/S 
v.s. H/Hs, while Figure |H1 shows the field dependence 
of M/S = {Mrj - SM)/S and SM/M^i, aU for S = i 
case. In both Figures the results for D — and D = 0.04 
are shown. 

The overall field-dependence of the considered quanti- 
ties is the same: quantum corrections are suppressed by 
the fieldi^ For the case of the pure Heisenberg model 
{D = 0) the ground state above the saturation field 
(H > Hs) is classical and all quantum corrections vanish 
ai H = Hg. Since in H _L D configuration the DM in- 
teraction prevents the full spin polarization by the field, 
quantum effects survive in the region H > H^. How- 
ever, with the further increase of the field they diminish 
quickly, see Figs. |Z|and|S| 

As is shown in Figures and |S1 different quantities are 
affected differently by the DM interaction. While we will 
elaborate on the detailed features of such effects below, 
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we would like to note that the quantum corrections to all 
the considered quantities are undoubtedly suppressed by 
the DM interaction compared to the D — case for the 
fields H <Hs. 

Besides the proliferation of the quantum fluctuations 
into the H > Hg region and an overall suppression of 
them for H < there are two notable features shown in 
Figs. [TjandlHl First, while changes in the l/^-corrections 
to the energy due to a modest D = 0.04 are hardly notice- 
able, changes in the quantum fluctuations of the on-site 
magnetization SS are substantial. Moreover, the rela- 
tive role of such DM-induced changes seems to be signif- 
icantly amplified by the applied field. For instance, at 
H = Hs/V^ SS is suppressed by almost 40% compared 
to its D = value due to the DM-term of the magni- 
tude only 4% of J, see Fig. Second, the behavior of 
6M/Mci v.s. H is non-monotonic with significantly dif- 
ferent zero-field limit of this quantity for the D = and 
D ^ cases, see Fig. |H1 

Uniform magnetization, H ^ 0. We would like to 
discuss the zero-field behavior of 5M /Md first. The rea- 
son for the choice of such a quantity is simple: the ratio 
of the quantum 1/5* part, to the classical magne- 
tization, M = S sin (/?, shows explicitly the relative role 
of the quantum correction in the uniform magnetization. 
In addition, in the pure Heisenberg £) = case all the 
terms in the 1 / iS-expansion of the uniform magnetization 
vanish identically as, H ^ 0. Thus, at D = and H ^ 
6M/Mci is also related to the l/^-correction to the trans- 
verse susceptibility: X = Xci — Sx = {Mci — 5M)/H [recall 
that sm. if = H/Hs when D = 0]. Thus, using Eq. lf?7|) 
and after some algebra the zero-field value of 6M/Mci for 
D — Q case can be written as: 



lim lim = — 7k ^ / — 

H^OD^O Mai 2S ^ V I 



7k 0.55115 



7k 



25* 



(29) 



compare with the expression for 5x in Ref. l2a- This 
numerical value is marked by the upper arrow pointing 
towards the left y-axis in Fig. |S1 One can see, however, 
that the DM interaction changes this quantity drastically. 
Neither the classical nor the quantum part of the uniform 
magnetization vanish in zero field anymore because of the 
weak ferromagnetism induced by the DM term. In other 
words, canting angle if does not go to zero at ^ 0, 
but approaches a finite value ~ D/2. This leads to the 
following expression for the opposite order of limits: 



lim lim — — 

D^O H^O Mr] 



1 



which is nothing but the expression for the 1/5- 
correction to the on-site (staggered) magnetization 53/ S 
in the pure Heisenberg model at if = 0. This value is 
marked by the lower arrow pointing towards the left y- 
axis m Fig. ISl Clearly, the behavior of SM /M^i v.s. H 
and D is singular as the limits H ^ and D taken 
in different order lead to different results. For D 



there is a maximum in SM/Md v.s. H around H ~ D/2 
and then this quantity decreases roughly parallel to its 
D = Q counterpart, see Fig. |H1 

On-site magnetization. To address the question on 
the sensitivity of the quantum component of the on-site 
magnetization, 5S, to the DM interaction it is useful to 
analyze the quantity AS^^ = SSq — SS as a function of 
D, where the subscript refers to the D = value. Using 
Eq. (|25|l one can write: 



AS 



DM 



SSn - SS 



S 



-y 



25^ v., 



Al A. 



I .k" 



where A-^ = C\ ^ C27k with C\ and C2 from Eq. l(T7|) 
and fk = 0^1^/45 from Eq. H18() as before, and and 
Al. arc their £> = hmits, Eq. lO. O ne can show that 
in the limit Z) ^ 1 the difference SSq — SS is dominated 
by the contributions from the region in k-space close to 
the (tt, 7r)-point with the rest of the Brillouin zone con- 
tributing to the subleading terms only. It may thus be 
concluded that the behavior of the DM-induced change 
in SS should be related to the DM-induced (tt, 7r)-gap dis- 
cussed in Sec. IIIBI 2. After some algebra one can obtain 
the asymptotic expressions: 



H = 



AS 



DM 



S 



2ttS 



^Dsm2ip, D-^H <H, , (32) 



where sin 1^9 ~ H/Hs. Thus, aside from a somewhat more 
subtle numerical coefficient for the result at H = Hs, 
the considered quantity follows almost identically the D- 
dcpendcnce of the (tt,7t) gap, compare with Eqs. H22|l . 
(|23|l . This asymptotic consideration also shows that ex- 
cept for the close vicinity oi H = and H = Hg the DM- 
induced change in SS is proportional to oc D^l^ . That is, 
it is significantly amplified compared to the other quan- 
tities and is not small already for the moderate values of 
D as we observed in Fig. [7| To demonstrate the validity 
of this asymptotic result, Eq. H32() . one can take the in- 
tegral in Eq. H31|) numerically and compare it with the 
results of Eq. 132|) in the D \ limit. Such a compar- 
ison for a representative value of the field H = Hs/\/2 
(sin2(p = 1 in Eq. (|32|l ) is shown in Fig. EJa) in both 
log-log and linear scales. The coefficient B = 0.48 for the 
subleading term {0{D)) is used in this plot. One can see 
an excellent agreement of the numerical integration of an 
exact expression (squares) with the asymptotic formula, 
Eq. (EH, (solid lines). 

Uniform magnetization. A similar analysis can be 
given to the other quantities. We have already discussed 
above the DM-induced changes in SM/Mci in the zero- 
field limit. Introducing AM^*^ = JMo - SM to charac- 
terize the effect of the DM interaction on the quantum 
component of the uniform magnetization, where SMq is 
the D = limit of SM , after some algebra one obtains 




D/J D/J 

FIG. 9: (Color online) The DM-induced changes in the l/S-corrections to (a) the on-site magnetization at H = Hb/V^ and (b) 
the uniform magnetization ai H — Ha- Exact results obtained by using the numerical integration in (a) Eq. I|25^ and (b) Eq. 
(I27I I (squares) and the asymptotic expressions (a) Eq. (EU and (b) Eq. ^ (sohd hnes) are shown, (a) A — l/27rS', B = 0.48, 
(b) A = 1/2,-rS. S = 1/2. 



the asymptotic expressions: 



AM 



DM 



1 

25^ 



r 0.1574, 

AiD, 



F = , 



(33) 



Clearly, these are less trivial dependencies than for the 
on-site magnetization, Eq. H32|) . The number for the 
H = case is simply the difference of the results in Eqs. 
(|29|l and (|30() discussed before. For the intermediate-field 
regime D <^ H < Hg the DM-induced changes in the uni- 
form magnetization are not dominated by any particular 
region in the k-space, but are accumulated over the en- 
tire BZ. Therefore, there is no straightforward relation 
between the D- or i?-dependence of the (tt, 7r)-gap and 
AM^*^. The expression for Ai{H) is rather cumber- 
some and we only list its value for a representative field: 
Ai{Hs/^/2) = 0.1580. On the contrary, in the regime 
H = Hs the DM-induced changes in SM/Md are domi- 
nated by the vicinity of the (tt, 7r)-gap, but only with the 
logarithmic accuracy, see Eq. since the subleading 

term is 0{D^/^). To demonstrate the validity of the lat- 
ter asymptotic expression we evaluated {6Mq — SM)/Mci 
numerically using Eq. (|27|) and compared it with the re- 
sults in Eq. H33() for H = Hg. An excellent agreement 
between the numerical (squares) and asymptotic expres- 
sions (solid lines) is shown in Fig. ^h). 

Energy. A similar consideration can be given to the 
l/S'-correction to the energy. For the intermediate-field 
regime D H < Hs one finds that {SEq - 5E) /iS = 
A3D + 0(1)3/2), gg^j^g SM^^' but with much smaller 
numerical coefficient. For example, for H = Hg/ the 



expression for the coefficient ^3 can be simplified to 

^3 4E(^-i)-0-0259. (34) 

At the saturation field the energy correction can be ap- 
proximately written as: 

which contains an extra power of D^^^ in comparison 
with SM/Mai, Eq. (Pt. 

Summary. Altogether, the quantum 1/S corrections 
are changed substantially under the influence of the DM 
interaction. Many of the considered changes depend on D 
non-analytically, see Eqs. (ESJ, and (|^ . and some 

exhibit a singular behavior as for the case of 6M/Mci in 
{H,D) -> limit, see Eqs. and Fig. El The 

most substantial change is seen in the on-site magnetiza- 
tion where for the broad field range SSq — SS(x D^^^, see 
Eq. (|32|l and Fig. [7\ Except for the vicinity of the satu- 
ration field where the quantum effects arc enhanced (or, 
rather, induced for H > H^), the quantum fluctuations 
are suppressed by the DM interaction, see Figs. [71and|Hl 



III. K2V3O8 
A. Experiments and Hamiltonian 

K2V3O8 has been studied experimentally by ma gne- 
tometry and neutron scattering, as reported in Ref. 122 . 
This material contains weakly coupled square lattice 
planes of spins S = The spins interact antiferromag- 
netically with the superexchange constant J = 12. 6K. 
In the ground state spins are pointing along the z-axis 
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while in a low external field {H < IT) two types of tran- 
sitions were seen. For the field along the z-axis the well- 
known "spin- flop" -type transition occurs at Hsf'- spins 
suddenly "flop" in the x — y plane and cant towards the 
field upon further increase of the field in a configuration 
similar to Fig. ^a). For the field H _L z an unusual 
"spin-rotation" transition is seen: spins perforin a rota- 
tion by 7r/2 angle gradually from along the z-axis into the 
x — y plane within a very small field range < H < Hsr- 
Above the "spin-rotation" field Hsr spins are in the x — y 
plane and cant towards the field in a configuration sim- 
ilar to Fig. n^b). Such a behavior made the authors of 
Ref. to conclude that there are small anisotropics of 
two types which accompany the superexchange in this 
material: the easy-axis anisotropy along the z-axis and 
the DM-anisotropy also directed along the z-axis. There- 
fore, K2V3O8 is described by our Hamiltonian (QJ with 
an additional easy-axis term: 



5nea=EY,S!S'- . (36) 



As is shown in Ref. "25, with an appropriate choice of pa- 
rameters the Hamiltonian ((T))-I- H36|l describes excellently 
all the observed transitions. 

As is mentioned in Sec. ^ although the easy-axis con- 
stant E is of the order of 0{D^), it competes with the 
DM-term on equal footing. While for an ideal staggered 
configuration of the DM vectors it is expected that the 
easy- axis constant is E = and the system possesses 

a hidden 0(3) symmetry^S^ a non-staggered component 
of the DM vector generally breaks such a high symme- 
try and leads to E < £)^/2i2i In realistic systems, like 
K2V3O8 the easy-axis term is, de facto, E > due 
to some other microscopic processes. We thus consider 
D and E as independent parameters, but still E ~ D^. 
In the rest of the Section we consider the experimentally 
relevant case E > /2 and discuss the case E < /2 
in the end. 

One can show that for the case of vector D being along 
the "easy" axis and for E > (we continue to use 

J = 1 units here), the easy-axis anisotropy overcomes the 
DM-term in the ground state. That is, instead of stay- 
ing in the plane _L D and canting towards each other, 
as favored by the DM interaction, spins align along the 
"easy" axis as dictated by the easy-axis term and as if 
the DM-term was absent. Therefore, in the ground state 
of this Hamiltonian a,t H — and E > D^/2 there is 
no weak ferromagnetism and the DM interaction is "hid- 
den". However, the external field helps to reveal the 
presence of the DM interaction in the system. It is the 
interplay of the easy-axis and the DM-term which causes 
an unusual spin rotation from the z-axis into the x — y 
plane in K2V30g. 



B. Ground state and excitation spectrum in a field 

We thus consider the Hamiltonian (fTT)+ H3fci|) as a start- 
ing point and take the ex per imental estimate for the DM 
constant D = 0.04, Ref. |22. For the easy-axis constant 
we use the value E = 0.0017, somewhat higher than es- 
timated in Ref. [2^ for the reasons discussed below. 

Before the technical discussion we would like to de- 
scribe the field evolution of the ground state from the 
symmetry point of view. In zero field for H < HspiHsR) 
the continuous 0(3) symmetry is broken by the easy-axis 
term to a discrete Ising-like one. Thus, all the excitations 
in such a phase should be gapped. For H|jz above the 
spin- flop field, H > Hsf, the continuous symmetry is re- 
stored to 0(2). This symmetry corresponds to the free- 
dom of choice of the direction of the spins' component in 
the x — y plane. Since the DM vector is also |jz, its pres- 
ence does not change this consideration except that the 
spins acquire a small canting angle towards each other 
in the x — y plane, see Fig. Q^a). Therefore, one of the 
magnon branches should be gapless in this phase. This 
is identical to the H||D field configuration discussed in 
Sees. Ill Al l. iriBI l without the easy-axis anisotropy. For 
the in-plane field H _L D||z above the spin- rotation field, 
H > Hsr, all the continuous symmetries remain broken 
and the ground state is identical to the one considered 
in Sec. Ill AI 2 and Fig. ^b) where spins are in the x — y 
plane and are canted towards the field. In such a phase 
no gapless excitation should exist. 

With these ideas in mind we will analyze classical spin 
configurations and excitation spectrum of the Hamilto- 
nian lfT|l + (|36|l . Because of the easy-axis term there are 
two different phases for each of the field directions: be- 
low and above the spin-flop (spin-rotation) transition for 
H||z (H _L z). 



1. Out-of-plane field. 

For the H|jD|jz configuration the classical energy of 
the Hamiltonian ((TT)+ 136|) can be written as: 

^cos9ACOseBil + E) 

— sin 0^ sin (cos 2(p -f I? sin 2(p) (37) 

- {COS0A + 00363) , 

where 9a and 6b are the angles with the z-axis made 
by spins in the A and B sublattice, respectively, ip is 
the angle in the x — y plane as in Fig. ^a). Using the 
symmetry arguments one can show that for the case E > 
D^/2 two choices of the spin configuration are possible: 
(i) 6a = 0, 6b = TT, and (u) 9a = 6b = 6, as in Fig. 
^a). Needless to say, these two choices correspond to 
H < Hsf and H > Hsf, respectively. For the first 
choice the spins are pointing along the z-axis and the 
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FIG. 10: (Color online) The field-evolution of the gaps for 
the spin-flop field configuration H||D||2. The spin-flop tran- 
sition occurs at Hsf, Eq. 14111 . Below the transition the gaps 
depend linearly on H, Eq. 14311 . Above the transition one 
of the gaps is zero (dashed) and the other corresponds to a 
uniform precession of the field-induced magnetization (solid), 
Eq. 14511 . The uniform precession mode for the E = D = 
case is shown by the dotted line (A = H). Inset: same results 
in the region of the small fields. Saturation field for K2V3O8 
is H, = 8SJ = 37.6T. 



classical energy is field- and DM-independent: 



E 



(1) 



2NS^ 



-{1 + E) 



(38) 



For the second choice the energy is minimized by the 
following choice of angles: 

1 U 2 

in = - tan-^ D ; cosO ^ — , (39) 

^2 8S 1 + E + VTTD^ 

where (p gives a small field-independent canting of spins 
in the x — y plane. For small fields spins are almost in 
the X — y plane as « 7r/2, all in a very close similarity 
with our discussion of the DM-only case, Sec. Ill Al 1 . 
Fig. n^a), and Eq. Q. Using these angles H39() one can 
rewrite the classical energy as: 



E 



(2) 



27V52 



From the condition E 



(1) 
•1 



lQS^-{l + E + y/\ + D^) 



(2) 

E\ we can now find the field 



• (40) 



at which the spin-fiop transition between the two spin 
configurations occurs. This gives: 



HsF^ 4.S\/2E -D'^ + E^ ^ AS\/2E - , (41) 

which is well-defined only if £^ > /2. Because E ~ 
<C 1 one can neglect E"^ , D^, etc. as being extremely 
small. One can verify that the first derivative of the 
energy with respect to field has a jump as is expected 



for the first-order transition. The angles also experience 
discontinuous jumps to finite values at Hsf as can be 
inferred from Eq. I|39|l . 

The field evolution of the excitation spectrum within 
the spin-flop problem has been studied in the past and is 
well-known (see, e.g., Ref. l3^ . In our case there are some 
quantitative changes due to the DM-term, but the overall 
behavior is very similar. Below the transition there are 
two branches: 



u>t=^S^{l + EY-^l{l + D^)±H , (42) 
with the gaps: 

I^^=HsF±H , (43) 

where the spin-flop field is given in Eq. H41|l. 

Above the spin-flop transition the spectrum can be 
written as: 



c^k = 45v/l + i?V(l + 7k)(l-<^i7k) , (44) 
where ..^ri + i. ^^^/^^^^ ^ ' 



1 + E + VTTd^J VTTd^ ' 

with one mode gapless, Atttt = 0, and the other mode 
having a gap which can be written as: 



1 + s + VTT^ ^ 



(45) 



This mode reaches the uniform precession behavior A ~ 
H Sit H ^ Hsf- All these field-dependencies of the 
gaps are shown in our Fig. 1101 One can see a jump in 
one of the gaps associated with the spin-fiop transition. 
The dotted fine is A = H {E = D = result). As we 
discussed above, one of the modes is necessarily gapless 
for H > Hsf, same as for the E = case. Sec. Ill Bl l. 



2. In-plane field. 

For the H _L D||z field configuration the spin arrange- 
ment within the spin-rotation phase < H < Hsb. is 
given in Fig. ^2 Note that ip is the angle between the 
projections of spins on the x — y plane and the x-axis, 
while 9 is the angle between the actual spin direction 
and the z-axis. One can see that the angle with the z- 
axis 9b — i^ — 9a for the spins in different sublattices. 
Both angles, (p and ^, are coupled non-trivially in the 
classical energy of the Hamiltonian ((T|)+ l|36() which can 
be written as: 



Ed 
2NS^ 



= -cos^9{l + E) 



(46) 



H 



— sin 9 (cos 2ip + D sin 2ip) sin 9 sin tp . 

Minimizing the energy, after some algebra, leads to the 
following expressions for the angles: 



tan (/J : 



E 
D 



sin 9 



H VD^ + 

AS 2E-D^ + E^ ' 



(47) 
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FIG. 11: (Color online) The projections of the spin arrange- 
ment on the X — Z-, X — y-, and y — z-planes within the spin- 
rotation phase < H < Hsr- While ip is the angle between 
the projections of spins on the x — y plane and the x-axis, 
is the angle between the actual spin direction and the z-axis. 
H is chosen along the y-axis. 



Somewhat surprisingly, the ip angle remains field- 
independent throughout the spin-rotation phase. Since 
the easy-axis interaction E is of the same order as ^ 
1, iy9 « E/D is also small. On the contrary, 9 changes 
drastically (but continuously) from at 7J = to 7r/2 at 
the critical field Hsr- The latter can be found from Eq. 
(HTjl as: 



H; 



SR 



AS 



2E- 



D 



(48) 



where we omitted 0{E'^) terms. 

For the fields above the transition, H > Hsr, a differ- 
ent solution minimizes the classical energy in Eq. (|46|1 : 
9 = 7r/2 is constant and f is field-dependent. That is, the 
spins are lying in the x — y plane as depicted in Fig. H^b). 
It is easy to see that ioi 9 — tt/2 Eq. 14t)|) reduces to the 
classical energy for the DM-only case, Eq. Q). Thus, the 
easy-axis term is effectively "switched-off" by the field for 
H > Hsr- Because of that the field-dependence of the 
canting angle (p also remains the same as in Sec. Ill AI 2. 
Eq. ©. 

Using the canting angles, Eqs. (|T7|l and ©, one can 
simplify the expressions for the classical energies below 
and above the spin-rotation transition to: 



2NS^ 
2NS^ 



-{1 + E) 



E 



(49) 



-1 



-(- 



JP 

1652 2E~D^ ' 

2 

D I , (50) 



where 0{E'^) terms are neglected. From these formulae 
one can verify that the spin-rotation transition is of the 
second order as the energy and its first derivative are con- 
tinuous through the transition. The angles also change 
continuously at Hsr as can be seen, for instance, for 
E,D^ <$:l: (p< ^E/D matches with (p> = D/2 + H/8S 
a.tH = Hsr, Eq. 



These canting angles can be translated directly into 
the field-induced uniform magnetization which has been 
measured as a function of applied field, see Ref. I^^. In 
small field the uniform magnetization depends linearly 
on field, but has different slopes for H < Hsr and H > 
Hsr- On the other hand, such slopes can be extracted 
from our Eqs. H49|) and (|50|l which give: 



dH 



fdAr 



V dH 



1 



D2 
2E 



(51) 



This relation is useful to extract the ratio /E because 
the slope in the magnetization v.s. field can be measured 
rather reliably and is the subject of less experimental 
error than, say, the transition field itself. It is using this 
relation and D = 0.04 we have extracted the value E — 
0.0017 for K2V3O8 used throughout this Section. 

After the easy-axis and the DM constants are chosen 
one can find the T = values of the spin-flop and spin- 
rotation fields from Eqs. (|41l) and (|48l) . Using the values 
of gc = 1.922 and gab = 1-972 for K2V3O8 reported in 
Ref. Hi one finds: Hsf « 0.85T and Hsr « 0.89T. 
These theoretical values are within the same range as 
H^g^ = 0.85T and TJ^a = 0-65T found in Ref. [H from 
the neutron scattering and magnetization measurements 
and H^g^ = 0.95T and Hg^^ = 0.55T estimated from 
the upturn in the thermal conductivity v.s. field in Ref. 
Is^ Since all these experimental results were obtained 
at a finite temperature of the order of the Neel order- 
ing temperature for this material, the differences between 
the theoretical results and the data are to be expected. 
The closer agreement with the data for the spin-flop field 
are due to the first-order nature of that transition. It 
is knownS that the critical field for the spin-flop tran- 
sition is only weakly temperature-dependent. For the 
spin-rotation transition, on the other hand, one can ex- 
pect a suppression of the critical field with the tempera- 
ture. Overall, the theoretical results for the critical fields 
are in a close accord with the experimental data. 

One can now address the field-dependence of the exci- 
tation spectrum for the spin-rotation problem. As is ex- 
plained in Sec. IIIBI we need the field-dependence of the 
canting angles, a transformation from the local to the lab- 
oratory frame for spins, and the subsequent quantization 
of the spin Hamiltonian to obtain the spectrum. While 
above the transition, H > Hsr, the first two steps are 
essentially identical to the ones performed for the DM- 
only case in Sec. IIIBI 2. below the transition, H < Hsr, 
this procedure requires an additional and a somewhat 
cumbersome algebra. We thus simply provide here the 
results for the gaps in the spin-rotation phase H < Hsr'. 



Aoo = W ^SDHsR { 1 + ^^M^] , (52) 



A^^^ JASDHsrII- 



^SR 



(53) 



where Hsr. is from Eq. (|48|l and we neglect 0{E'^) terms 




FIG. 12: (Color online) Aoo/4S'J v.s. H/Hs for the in-plane 
field direction, D = 0.004, and two values oi E: E ^ 0.0017 
(solid), Eqs. ^ and and £ = (dashed), Eq. (HHJ. In- 
set: same plot in the region of small fields. The spin-rotation 
transition is marked by the dotted vertical line. Saturation 
field for K2V3O8 is Hs = 8SJ = 37.6T. 



FIG. 13: (Color online) A^^/ASJ v.s. H/Hs for the in-plane 
field direction, D = 0.004, and two values oi E: E = 0.0017 
(solid), Eqs. (|H and (EHJ, and £ = (dashed), Eq. J^TJ- 
Inset: same plot in the region of small fields. Saturation field 
for K2V3O8 is Hs = 8SJ = 37.6T. 



as before. 

For the fields above the transition, H > Hsr, the 
changes with respect to the DM-only case. Sec. IIIBI 2. 
are rather incremental and we can afford listing the pa- 
rameters of the spin- wave Hamihonian: 



Ci = 1 + D tan ip , 
C2 + C3 = — cos 2(p — D sin 2tp , 
- C3 - 1 + S , 



(54) 



which one can compare to an equivalent expression in Eq. 
(|17|l . Thus, the easy-axis term leads to rather obvious 
changes in the spin- wave energy, Eq. H18|) . The gaps are 
changed accordingly and are given by: 



Our results for the field-dependence of the Aqo and 
Atttt gaps, both below and above the spin-rotation tran- 
sition, together with the results for the DM-only {E = 0) 
case are shown in Figures El and El respectively. One 
can see that for the Aqo gap above the transition the 
results for K2V3O8 are essentially indistinguishable from 
the DM-only case. It can be shown that the difference be- 
tween them scales as esc 0{EH) = 0{D-^H). The changes 
in Atttt with respect to the DM-only case are more sub- 
stantial. They scale as 0{E/D) for the fields H > Hsr 
but diminish closer to the saturation field. Thus, the 
major effect of the easy-axis anisotropy on the excitation 
spectrum is the effective "off-set" of the external field by 
Hsr + 4:3 D, as can be seen from comparison of Eqs. l|57|l 
and 



Aoo^AsJ^sinipn+^DtauLp + E] , (55) 



A^^ = 4:Sy/2{Dt£intp- E) 



(56) 



X J n + Z) tan - sin ^ ) , 



which should be compared to the DM-only results, Eqs. 
(|19|l and H21|l . respectively. It is important to note that 
in the vicinity of the transition field the (tt, 7r)-gap has 
the following non-analytic field-dependence: 



A,, = ^ASD{H - Hsr) 



(57) 



for H > Hsr- This is reminiscent of Eq. H22|l for the 
H ^ Hs case, but with the field "shifted down" by Hsr+ 
ASD. 



E < 12 case. 



For completeness, we briefiy discuss here the case E < 
In such a case the ground state and the spectrum 
field-dependencies are qualitatively the same as in the 
DM-only case ^ considered in Sees. Ill Al and III Bl In 
fact, for the most interesting H _L D configuration the 
canting angle in the ground state is simply identical to 
Eq. (jHJ because the easy-axis anisotropy is effectively 
switched-off. In particular, the zero- field ground state 
for E < D^/2 possesses the same uniform magnetization 
as if the easy-axis term would be absent. Interestingly 
enough, the high- and intermediate-field behavior of the 
(tt, 7r)-gap also depends only on the DM-term, Atttt cx; 
{DHf/'^ for H :^ D. On the other hand, the zero-field 
gap in this case is reduced: At^-k = 43 V — 2-E. Thus, 
the zero-field canting angle ip D/2 is not in one-to-one 
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correspondence anymore with the zero-field (tt, 7r)-gap, in 
contrast to the DM-only case, model QJ. 

A particular case E — is also interesting as it 

shows no gap in zero-field, yet a finite-field behavior of 
the gap is governed by the DM-term: A^tt oc {DHY^^, 
the result obtained in Refs. Another interesting 

feature of that case is that the classical energy in zero- 
field, Eq. (|46|) . is degenerate for any choice of 6, while 

= ^ tan^^ D ~ D/2. That is, the weak ferromagnetism 
with the uniform magnetization of an arbitrary strength 
from the range < AI < SD/2 can occur in zero- field 
ground state as a result of the spontaneous choice of 6 in 
the spin configuration depicted in Fig. 1111 This scenario 
has been also discussed on Ref. l25l 

4- Summary. 

Altogether, an additional easy-axis anisotropy intro- 
duces some new transitions in the physics of a realistic 
system in comparison with the DM-only case, but leaves 
most of the results for the field above the corresponding 
transition unchanged. That is, for the out-of-plane field 
configuration one of the gaps in the spin-wave spectrum 
behaves similarly to a uniform precession mode, A « _ff, 
while the other mode stays gapless. For the in-plane field 
both modes are gapped. One of them is the uniform pre- 
cession mode while the other depends non-analytically on 
the field and shows a non-trivial behavior in high field. 
We would like to note that some indications of such a 
drastically different behavior of the gaps for different field 
directions have been observed in the field-dependence of 
the thermal conductivity'^'^ in K2V3O8, see Ref. |3 One 
of the crucial differences of the K2V3O8 spectrum field- 
dependence from the idealized DM-only case, Eq. 
is the effective off-set of the field for the in -plane field 
direction, leading to the oc D{H — Hsb.) behavior of 
the gap. Thus, in a realistic system, the low-field behav- 
ior of th e DM-induced gap is more delicate than simple 
cx y/DH suggested recentlyiSiiiiiS 

IV. CONCLUSIONS 

We conclude by summarizing our results. We have 
studied the effects of external field on the properties of an 



ordered Heisenberg antiferromagnet with the DM inter- 
action. Utilizing the spin- wave theory we have calculated 
the quantum correction to the energy, on-site magnetiza- 
tion, and imiform magnetization v.s. H and D. We have 
shown that the spin-wave excitations exhibit an unusual 
field-evolution of the gaps which leads to various non- 
analytic dependencies of the quantum corrections. We 
have found that the gap in one of the magnon branches 
follows a particularly interesting behavior v.s. field for a 
specific field direction: A^^r is oc D, esc D^/^, or oc D^/'^ 
depending on the field. A similar behavior is reflected 
in some of the quantum corrections as well. We have 
demonstrated that for the fields < H < Hg the DM in- 
teraction suppresses quantum fluctuations. The on-site 
ordered moment is affected most strongly by such a sup- 
pression. It is only in the regime of the field close to 
the saturation where the DM term enhances quantum 
fiuctuations, effectively leading to a proliferation of the 
quantum effects into the classical saturated phase. We 
have shown that some quantum corrections are singular 
in the limit of H,D ^0. 

We also have considered K2V3O8 which is well de- 
scribed by the spin-i, Heisenberg AF on the two- 
dimensional square lattice, with the additional DM and 
easy- axis anisotropics. We have demonstrated that 
K2V3O8 is an excellent candidate for the observation of 
the unusual field-dependence of the gaps, characteristic 
to the 2D and 3D AFs. In fact, an additional easy-plane 
anisotropy mak es it potentiall y possible to observe a non- 
analytic A oc ^JD{H — Hsb.) behavior of the gap. 
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